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Photo-disintegration cross section measurements on 186 W, 187 
implications to the Re-Os cosmochronology 
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Cross sections of the 186 W, 187 Re, 188 Os(7, n) reactions were measured using quasi-monochromatic 
photon beams from laser Compton scattering (LCS) with average energies from 7.3 to 10.9 MeV. 
The results are compared with the predictions of Hauser-Feshbach statistical calculations using 
four different sets of input parameters. In addition, the inverse neutron capture cross sections 
were evaluated by constraining the model parameters, especially the El strength function, on the 
basis of the experimental data. The present experiment helps to further constrain the correction 
factor F„ for the neutron capture on the 9.75 keV state in 187 Os. Implications of F a to the Re-Os 
cosmochronology are discussed with a focus on the uncertainty in the estimate of the age of the 
Galaxy. 
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I. INTRODUCTION 



Owing to the long half-life of 187 Re, the 187 Re- 187 0s 
pair may serve as a cosmochronometer to measure the 
duration of stellar nucleosynthesis that precedes the so- 
lidification of the solar system p| . By adding the age of 
the solar system (~4.6 Gyr), it provides the age of the 
Galaxy. The facts that both 186 Os and 187 Os are pro- 
duced only by the s-process nucleosynthesis apart from 
the cosmoradiogenic yield of 187 0s and that the isotopic 
solar abundance ratio of 186 Os and 187 Os 0] is available 
make this chronometer potentially reliable in the sense 
that it is independent of r-process models. The quantita- 
tive interpretation is, however, complicated by the pos- 
sible enhancement of 187 Re- 187 0s transmutation rates in 
stellar condition, the stellar production and destruction 
of 187 Re and 187 0s during the chemical evolution of the 
Galaxy, the possible existence of s-process branchings at 
185 W and 186 Re, and the neutron capture by the 9.75 
keV first excited state in 187 0s US 110. 



The last issue on the effect of neutron capture on the 
9.75 keV state in 187 0s was raised in Ref. 0. In the 
local approximation, the ratio of the s-process yields of 
186 0s and 187 0s, N s ( 186 0s)/N s ( 187 0s), can be expressed 



as 



7V s ( 187 0s) 



F r 



(a)( 186 0s) 



N s ( 1S6 Os) <T (cr)( 187 0s) 



(1) 



Here, (cr)( 186 0s) and (er)( 187 0s) are the Maxwellian- 
averaged neutron capture cross sections on 186 Os and 
187 0s in the ground states, respectively. The F a value 
accounts for the correction to the cross section due to the 
neutron capture on the 9.75 keV state in 187 0s which is 
substantially populated at typical s-process temperatures 
T -1-3 x 10 s K. It is defined by 



F„ = 



(q)( 187 0s) 
(ct)*( 187 0s) 



(2) 



with (ct)*( 187 0s) being the Maxwellian-averaged neutron 
capture cross section on 187 0s at a given stellar tempera- 
ture. Here, the first excited state in 186 Os and the second 
excited state in 187 0s lie at higher excitation energies of 
137 keV and 74 keV, respectively, and therefore their 
contributions to neutron capture in the stellar condition 
may safely be ignored. In 1970s, it was of critical con- 
cern whether or not the F a value exceeds unity, because 
it has a great impact on the age of the Galaxy; the larger 
the F a , the smaller the age. However, there was a large 
spread in the early estimate (0.80-1.10 0, 0.8 Q, and 
-1.5 0). 



2 



This concern, combined with the fact that a direct 
measurement of neutron capture on the 9.75 keV state is 
virtually impossible, led to measurements of neutron cap- 
ture on 186 > 187 < 188 Os 0,0] and neutron inelastic scatter- 
ing to the 9.75 keV state in 187 0s [H[l3- The measured 
neutron capture cross sections were in good agreement 
with those of the earlier measurements [3, UM- But, 
the statistical analysis of the capture data gave a lower 
bound of 0.30 b to the inelastic scattering cross section 
a nn i at a neutron energy of 30 keV, and by combining 
with an upper limit of a nn i = 0.5 b 15], gave F a ~ 1. 
In contrast, the two neutron inelastic scattering data, 
a nn , = 1.13 ± 0.2 b at 60 keV and 1.5 ± 0.2 b at 
34 keV 0, are consistent with F a =0.80-0.83 and 0.80, 
respectively, within the statistical models. 

More recent efforts have been made toward a unified 
statistical model analysis of all available data includ- 
ing measurements of (n, 7) cross sections |l6j and elas- 
tic/inelastic scattering cross sections 0] on a neighbor- 
ing nucleus 189 Os in which the ground state with J 71 ' = 
3/2~ and the first excited state with l/2~ at 36 keV ap- 
pear in the reverse order of the corresponding states in 
187 0s. These analyses showed that F CT =0.79-0.83. 

In the Hauser-Feshbach statistical model calculations, 
however, large uncertainties may arise from the 7-ray 
transmission coefficients rather than the neutron optical 
potential and the level density. Information on the 7-ray 
transmission coefficients for neutron captured states in 
the low-energy tail of the giant dipole resonance in 188 Os 
can be obtained in the inverse photo-disintegration of 
188 0s. But, the conventional Lorentzian model based on 
the previous photo-disintegration data on 188 Os may 
not be satisfactory for two reasons. One reason is that 
the data taken with 7-ray beams from the positron an- 
nihilation in flight exhibit non-vanishing cross sections 
even below the neutron threshold (see sec. IIIIBIl . The 
non- vanishing cross sections may be attributed to contri- 
butions from the positron bremsstrahlung. The other is 
that microscopic models can predict the El 7 strength 
function more reliably than the Lorentzian model. 

Besides F a , the effect of the s-process branchings at 
185 W and/or 186 Re was parameterized as Fb and inves- 
tigated within the framework of the schematic s-process 
models |4| . More recently, neutron capture cross sections 
were measured for neighboring 185 Re and 187 Re nuclei 
to derive the statistical model parameters from a con- 
sistent systematics [l9j. With the improved parameters 
for s-process analysis, a stellar model calculation for low- 
mass AGB stars showed that the local approximation was 
disturbed by the branchings at 185 W and 186 Re. How- 
ever, the precise physical conditions of the AGB model 
need to be scrutinized before any definite conclusion on 
Fb is drawn. This is clearly beyond the scope of the 
present study. Instead, the relevant statistical parame- 
ters, particularly the El 7 strength function, can further 
be improved by the photo-disintegration measurement on 
186 W. 

In the present study, we have measured photoneutron 
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FIG. 1: An observed spectrum of the LCS photons measured 
with an HP-Ge detector (solid line). The original LCS pho- 
tons simulated by the Monte Carlo code EGS4 is also shown 
with dashed line. The neutron separation energies S„ of 186 W, 
187 Re and 188 Os are indicated by solid lines. Fractional pho- 
tons with energies higher than S„ are responsible for the pho- 
toneutron reactions. 

cross sections of 186 W, 187 Re and 188 0s using tunable 
quasi-monochromatic 7-ray beams from laser Compton 
scattering (LCS). The photo-disintegration data allow us 
to constrain F a values within the Hauser-Feshbach sta- 
tistical model and discuss their implications to the Re-Os 
cosmochronology. Part of the present data has already 
been published |2Cj . 



II. EXPERIMENTAL PROCEDURE 

Photoneutron cross section measurements on 186 W, 
187 Re and 188 Os were performed at the National Institute 
of Advanced Industrial Science and Technology (AIST). 
Tunable quasi-monochromatic photon beams were gen- 
erated by Compton scattering of laser photons, with rel- 
ativistic electrons circulating in the storage ring TERAS 
plf. A Nd:YLF Q-switch laser at a wavelength of 527 
nm in second harmonics was operated at a frequency of 
2 kHz. The electron energy was varied in the range from 
450 to 588 MeV to produce LCS photons with the average 
energy from 7.3 and 10.9 MeV. A 20 cm lead collimator 
with a small hole of 2 mm in diameter was placed at ap- 
proximately 6 m downstream from the interaction area 
which defines a scattering cone of the LCS photons. The 
typical energy resolution was 10 % in FWHM. Further 
details on the experimental setup can be found in Ref. 

Figure ^ shows an energy spectrum of the LCS pho- 
tons measured by an HP-Ge detector with a relative ef- 
ficiency of 120 %. The energy was calibrated at 1460.8 
and 2614.5 keV with 40 K and 208 T1 radioactive isotopes 
of natural origin. A Monte Carlo simulation was per- 
formed with the EGS4 code |2^| to analyze the response 
of the Ge detector. The energy distribution of incident 
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FIG. 2: Pile-up spectrum of the LCS photons measured with 
a large volume (8"xl2") Nal(Tl) scintillation detector. In 
this spectrum the average photon number per laser pulse is 
5.4 which leads to the LCS photon intensity of 1.08 xlO 4 /s. 

LCS photons (dotted line in Fig. was determined so 
as to reproduce the observed response function (solid line 
in Fig. H). Photons with energies higher than the neu- 
tron threshold (5„=7.19 MeV for 186 W, S n =7.36 MeV 
for 187 Re and S„=7.99 MeV for 188 0s) are responsible 
for the (7,n) reactions. The fraction of these photons in 
the total photon flux and the average photon energy were 
obtained from the original LCS photon spectrum. 

The beam current of the electron storage ring decreases 
exponentially in a normal condition with a lifetime ~ 6 
hours. In the EGS4 Monte Carlo simulation, it was found 
that the electron beam size in the region of the interaction 
with laser photons varied with time: for example, from 
2.2 mm in diameter at 166 mA to 1.2 mm at 72 mA in the 
187 Re measurement. A space-charge effect is considered 
to be a main cause for the decrease in the beam size. 
This beam size effect, which was evident in long runs 
near the neutron thresholds, introduced uncertainties in 
the fraction of the LCS photon beam above the threshold. 
The resultant uncertainty was estimated to be 1 - 6 % in 
the present experiment. On the other hand, the average 
photon energy was determined well within 40 keV. 

The number of LCS photons was monitored during the 
experiment using a large volume (8"xl2") Nal(Tl) scin- 
tillation detector placed behind targets. A typical pile- 
up spectrum is shown in Fig. [5] The pulse height of the 
spectrum is proportional to the number of LCS photons 
per beam pulse. The photon flux was determined with 
3 % uncertainty based on a statistical analysis on the 
pile- up spectrum [24| . 

Metallic powders of 1246 mg 186 W, 996 mg 187 Re, and 
693 mg 188 0s enriched to 99.79, 99.52 % and 94.99 %, 
respectively, were pressed to self-supporting tablets with 
a diameter of 8 mm. The 188 Os powder included ma- 
jor contaminants of 189 0s (2.55 %), 190 Os (1.27 %) and 
192 Os (0.97 %). These tablets were mounted inside thin 
containers made of aluminum, and were irradiated with 
the LCS photon beams. The threshold energy of the 



(7, n) reaction on 27 A1 is 13.06 MeV, which is higher 
than those for 186 W, 187 Re and 188 0s. The present pho- 
toncutron cross section measurements were performed at 
energies below the threshold energy and therefore undis- 
turbed by the 27 A1(7, n) reaction. Further, measure- 
ments with an empty aluminum container (blank tar- 
get) showed that no background neutrons were produced 
from photo-disintegration of possible impurities in the 
aluminum. 



Emitted neutrons were detected by sixteen 3 He pro- 
portional counter (EURISYS MESURES 96NH45) em- 
bedded into a polyethylene moderator. Two sets of eight 
counters were placed in double concentric (inner and 
outer) rings at 7 and 10 cm from the beam axis. Time 
correlations (Fig. |3J) between the neutron signal and 
the laser pulse were measured to estimate the number 
of background neutrons that arrived randomly at the 
3 He detectors. These background neutrons were most 
likely produced by bremsstrahlung arising from collisions 
of electrons with residual gaseous molecules in the stor- 
age ring. In the moderation time distribution, constant 
events above 400 /xs and at small correlation times were 
taken to be background neutrons. The constant back- 
ground was further confirmed by using a 1 kHz laser and 
a wider (1 ms) time range (see, for example, Fig. 3 of 
Ref. H3). The background subtraction is included in the 
statistical uncertainties through the error propagation. 



The neutron detection efficiency was measured at the 
average neutron energy 2.14 MeV with a standard 252 Cf 
source. The dependence of the efficiency on neutron en- 
ergy was determined by a Monte Carlo MCNP simula- 
tion with the statistical accuracy less than 0.5 % [25| . 
The so-called ring ratio between the count rates of inner 
and outer rings was used to determine the average energy 
of emitted neutrons [23 • A polynomial fit to the energy 
dependence of the ring ratios was made as in , where 
an asymptotic value of the ring ratio (8.0) at 10 keV sim- 
ulated for the present neutron detector [25| was used as 
a constraint. The ring ratio varied between 7.4 and 2.9, 
indicating the neutron energies of several tens to a few 
hundreds keV for measurements close above the thresh- 
olds and up to 0.82 MeV at higher photon energies. The 
uncertainty in the neutron energy thus determined was 
estimated to be 10-15 keV, resulting in the uncertainty 
less than 0.8 % in the total neutron detection efficiency. 
Note that the total efficiency is nearly constant (44.4-44.3 
%) over the neutron energy from 1 MeV to 400 keV and 
that it slowly decreases to 39.5 % at 50 keV. The overall 
systematic uncertainty for cross sections was estimated 
to be 5.9-8.3 %, which was determined by the neutron 
emission rate of the calibration source (5 %), the number 
of the incident LCS photons (3 %), and the beam size 
effect. 
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FIG. 3: Typical time-correlated spectra between laser pulses 
and neutrons detected in the inner (left panel) and outer 
(right panel) rings. The time correlation was measured by 
using a TAC module with the neutron signal and the 2 kHz 
laser pulse being the start and stop, respectively. 



III. RESULTS 

A. Data reduction 

The photoneutron cross section measured with a 
monochromatic photon beam is given by 



a = 



N t N 7 fe n (E n ) 



(3) 



where n n is the number of neutrons detected with the 
3 He counters, Nt is the number of target nuclei per unit 
area, iV 7 is the number of incident photons, / is the cor- 
rection factor for a thick-target measurement, and e(E n ) 
is the neutron detection efficiency. The correction factor 
is given by / = (1 — e~^ d )/(nd) with the linear attenu- 
ation coefficient of photons, /i, and the target thickness, 
d. The attenuation coefficient was taken by interpolation 
from (2(| for the average energy of the LCS photons. The 
energy spread of the LCS photon beam in the full-width 
at half maximum makes negligible contributions (< 0.2 
%) to the determination of the correction factor /, which 
deviates from unity by no more than 6 % in the present 
measurements. 

Recently, a methodology was developed to deter- 
mine cross sections for reactions induced by a quasi- 
monochromatic photon beam [2(j. When the photon 
beam has an energy distribution of n^{E), N 7 a in Eq. @ 
has to be replaced by the integral / n 1 (E)a{E)dE . By 
expanding the cross section <r(E) in the Taylor series at 
the average photon energy Eq, the first term in the Tay- 
lor series a(E ) (the cross section at the average energy) 
was numerically evaluated along with the higher-order 
terms. The new methodology determines <j{Eq) in the 
energy region of astrophysical importance near threshold 
within 6 % corrections from the monochromatic approxi- 



mation (Eq. (3)). Photoneutron cross sections presented 
in this paper are based on this methodology. 



B. Photoneutron cross sections 

Photoneutron cross sections measured for 186 W, 187 Re 
and 188 Os are shown in Fig. Hand Table D The contri- 
butions from the reactions with the main contaminants 
( 189 0s, 190 Os and 192 0s) of the 188 0s target were esti- 
mated from the previous data from Ref. |27| . The error 
bars in Fig. include both the statistical and systematic 
uncertainties. Previously, cross section data were taken 
for 186 W and 188 Os [13] with quasi- monochromatic 
photons from positron annihilation in flight, as shown in 
Fig. 0] for comparison, fn addition, data for 186 W and 
187 Re were taken with bremsstrahlung [28L I2H l30j . In 
Ref. [3(j, yield curves obtained in small increments of 
the electron beam energy with 1-MeV spacing were con- 
verted to cross sections through an unfolding procedure 
based on the Penfold-Leiss method [3l]]. As mentioned 
in Ref. |3(j , it is a known fact [3^] that cross sections are 
not obtained correctly because of swings of the solution, 
though uncertainties resulting from the swings may be 
reduced significantly. 

The present 188 Os data are in reasonable agreement 
with the previous data [23] except at the low energy 
where the previous data exhibit non- vanishing cross sec- 
tions even below the neutron threshold of 7.99 MeV. The 
non-vanishing cross sections may be attributed to left- 
over in the subtraction of contributions of the positron 
bremsstrahlung admixed with the positron annihilation 
photons. 

The (7, n) cross section exhibits the threshold behavior 

mm of 



(4) 



Here, p is related to the neutron orbital angular momen- 
tum £ through p = i + 1/2. The values of p = 0.5 and 
1.5 are expected for the s- and p-wave neutron decays, re- 
spectively. The best fits to the present experimental data 
gave p = 0.47 and ctq = 78 mb for the 186 W(7, n) reac- 
tion, p = 0.67 and <tq = 117 mb for the 187 Re(7, n) reac- 
tion, and p — 0.53 and 00 = 143 mb for the 188 0s(7, n) 
reaction. The results show a rather pure s-wave char- 
acter for the 186 W and 188 0s(7,n) reactions, and sug- 
gest an admixture of p- and/or d-wave neutron emis- 
sions following the El excitation of 187 Re. The cross- 
section parametrization (Eq. was also made in the 
bremsstrahlung measurements [2$, |2j| ; the results are in 
rough agreement with our new experimental data. 
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FIG. 4: Photoneutron cross sections of the reac- 
tions 186 W( 7 ,n) 185 W (a), 187 Re( 7 ,n) 186 Re (b), and 
188 Os(7, n) 187 Os (c) extracted from the present data are plot- 
ted with filled circles. The previous data for 186 W [H and 
188 0s |27fl taken with quasi-monochromatic photon beams are 
also shown with open circles for comparison. Given with solid 
line is the best fits to the experimental data in Eq. 01 



IV. COMPARISON WITH THEORY 

A. Theoretical framework 

The cross sections measured in the present work 
are now compared with the predictions of the Hauser- 
Feshbach (HF) compound nucleus theory 0, HH- The 
uncertainties involved in HF cross section calculation are 
known not to be related to the theory of compound nu- 
cleus emission itself, but rather to the uncertainties asso- 
ciated with the evaluation of the nuclear properties en- 
tering the calculation of the transmission coefficients. It 
is therefore of prime importance to compare the effects 
of different nuclear inputs to estimate the reliability and 
accuracy of the predictions, especially when considering 
the reverse rates, i.e, in the present case, the radiative 
neutron capture rate which might be sensitive to differ- 
ent input parameters than the rate measured. In the 
present work, the nuclear level densities are derived from 
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-Eph is the average energy of the LCS photon beam. Cross 
sections are given as a ± Afr(stat.) ± A<r(syst.) in units of mb 
where Acr(stat.) and Acr(syst.) represent the statistical and 
systematic uncertainties. 

Nucleus 



3 W 



7 Re 



s Os 



ph 

(MeV) 


(J 

(mb) 


A<r(stat.) 
(mb) 


A<r(syst 
(mb) 


7.26 


7.40 
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0.43 


7.33 


15.7 
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1.0 
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14.3 


0.5 


0.9 


7.51 


19.6 


0.6 
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7.72 


23.7 


0.7 
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7.79 


21.1 


0.8 


1.3 


7.88 


25.0 


0.6 


1.5 


8.20 


27.6 


0.9 


1.6 


8.63 


35.7 


0.8 


2.1 


8.67 


35.5 


0.9 


2.1 


8.91 


38.8 


0.7 


2.3 


9.34 


46.3 


0.8 
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a The uncertainty includes those associated with the 
neutron detection efficiency (5 %), the photon flux (3 %) 
and the beam size effect (1 - 6 %) added in quadrature. 



two models, either the widely used back-shifted Fermi 
gas (BSFG) model based on the global parametrization 
of W or the microscopic calculations taking into ac- 
count the discrete structure of the single-particle spec- 
tra associated with Hartree-Fock+BCS (HFBCS) poten- 
tials [37|] ■ This model has the advantage of treating 
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shell, pairing and deformation effects consistently, and 
for practical applications, has been renormalized on ex- 
isting experimental information (low-lying levels and s- 
wave neutron resonance spacings whenever available as in 
the cases considered here). The transmission coefficients 
for particle emission is calculated either with the so-called 
JLMB semi-microscopic potential of j33] derived from the 
Briickner-Hartree-Fock approximation based on a Reid's 
hard core nucleon-nucleon interaction, or with the global 
phenomenological mass- and energy-dependent potential 
of Woods-Saxon type developed by [33 • 

The photon transmission function of particular interest 
in photoemission data is calculated assuming the domi- 
nance of dipole transitions in the photon channel. The 
electric- and magnetic-dipole (GDR) transition strength 
functions are usually described by a Lorentz-type func- 
tion where the energies and widths are determined by 
experimental data, whenever they exist, or by appropri- 
ate parametrizations. However, the calculation of the ra- 
diative capture or photoabsorption at low energies (and 
particularly in stellar conditions where the excited states 
of the target nucleus are thermally populated) is partic- 
ularly sensitive to the low-energy tail of the GDR of the 
compound system. The shape of the GDR is expressed 
most frequently by a generalized energy-dependent- width 
Lorentzian function adjusted on low-energy data [40| . To 
test such models, we consider here the Hybrid model ^l| 
which couples the GDR Lorentzian description at high 
energies with an analytical approximation to the theory 
of finite Fermi systems at energies below the neutron sep- 
aration energy |42|. 

In addition to the Hybrid model [ilf. the Quasi- 
Particle Random Phase Approximation (QRPA) model 
of 01 i s a l so considered here for estimating the pho- 
ton transmission coefficients. These QRPA calculations 
are self-consistently built on a ground state derived with 
the HFBCS approximation. The final -El-strength func- 
tions is obtained by folding the QRPA strength with a 
Lorentzian function to account for the damping of col- 
lective motions and the deformation effects. This global 
calculation based on the SLy4 Skyrme interaction has 
been shown to reproduce relatively well photoabsorp- 
tion and average resonance capture data at low energies 
pi H3. Note that both the Hybrid and QRPA mod- 
els differ not only in the predictions of the position and 
width of the GDR, but also in the energy dependence 
of its tail, which is an important quantity to derive the 
reaction rate. Since we are here mainly concerned with 
the GDR tail at low energies and the prediction of the 
reverse neutron capture cross section, both El-strength 
functions are renormalized on the available experimen- 
tal information on the position of the GDR peak and 
the corresponding maximum absorption cross section. In 
the case of the HFBCS+QRPA model, this adjustment 
is achieved within the folding procedure introduced to 
account for damping and deformation effects. 



B. Comparison between experimental and 
theoretical rates 

The final nuclear inputs considered in the present anal- 
ysis are summarized in Table [H] Four different sets are 
used to estimate the photoemission cross section as well 
as the reverse radiative neutron capture cross sections. 
The comparison between these 4 sets allow us to esti- 
mate the sensitivity of the cross sections to the various 
input quantities, but also the uncertainties affecting the 
final prediction of the neutron capture rate of astrophys- 
ical interest. In Fig. [5] our new experimental data are 
compared with the theoretical photoemission cross sec- 
tions. Also shown are previous measurements obtained 
in the vicinity of the GDR peak energy [2?], H3 • Most of 
the calculations agree relatively well with experimental 
data, down to energies close to the neutron threshold. 
The neutron optical potential and nuclear level densities 
influence the photoneutron cross section only in a small 
energy range of no more than 1 MeV above the neutron 
threshold, so that the global behavior of cross section is 
almost entirely dictated by the El-strength. Some spe- 
cific comments can be made for each reaction: 

• In the 186 W(7,n) 185 W case, the HFBCS+QRPA 
model predicts some extra strength at energies 
around 7.5-10 MeV with respect to the Hybrid 
model. This extra strength is clearly seen exper- 
imentally below 8 MeV but not above. However, 
all models overestimate the 7mb cross section at 
E = 7.26 MeV. 

• In the case of the 187 Re(7, n) 186 Re reaction, the 
low-energy data can only be reproduced when 
adopting the BSFG model of nuclear level densi- 
ties. The microscopic HFBCS-based model fails to 
describe the fast rise of the cross section at the neu- 
tron threshold. Around 11 MeV, the Hybrid and 
HFBCS+QRPA strength predict a relatively differ- 
ent cross section, the former one being compatible 
with the Goryachev et al. [3(| and the later with 
our more accurate measurements. 

• As far as 188 0s(7, n) 187 0s is concerned, all HF cal- 
culations reproduce relatively well the data, though 
the Hybrid model gives a lower cross section in the 
8.5-10 MeV energy range. 



C. Determination of the neutron capture cross 
sections 

We now estimate the stellar Maxwellian-averaged neu- 
tron capture cross section (a)* of astrophysics interest 
on the basis of the calculations presented above, i.e con- 
strained by the reverse photo-disintegration rate compat- 
ible with the new measurements. It should be recalled 
that the neutron capture cross section at energies of a 
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TABLE II: Overview of the 4 sets of nuclear ingredients adopted in the Hauser-Feshbach calculations. 



Label 


Level density 


7-strength 


Optical potential 


INP-1 


HFBCS [37] 


HFBCS+QRPA [43] 


JLMB [38] 


INP-2 


BSFG [36] 


HFBCS+QRPA [43] 


JLMB [38] 


INP-3 


BSFG [36] 


HFBCS+QRPA [43] 


Woods-Saxon [39] 


INP-4 


BSFG [36] 


Hybrid model[41] 


JLMB [22] 



Present 

Berman et al. (1979) 
-INP-1 



s 1 °- 10 



0.01 



,_, 0.10 




0.01 



,_, 0.10 



0.01 



■ Present 

• Goryachev et al. (1 973) 

INP-1 

INP-2 

INP-3 

INP-4 





■ Present 

• Berman et al. (1979) 

INP-1 

INP-2 

INP-3 

■-■-INP-4 



188 Os(y,n) 187 Os 

12 



9 10 11 

E [MeV] 



FIG. 5: Comparison between measured and calculated 
photoneutron cross sections for 186 W(7,n) (upper panel), 
187 Re(7,n) (middle panel) and 188 0s(7,n) (lower panel). The 
different theoretical predictions correspond to the input de- 
fined in Table El 



few tens of keV is mainly sensitive to photon transmis- 
sion coefficient at an energy close to and even below the 
neutron separation energy. For this reason, only the cal- 
culations that reproduce the photoemission rate at the 
neutron separation energy relatively well are retained. 



4 i 



3 - 



187 




Os(n,y) 188 Os 

Winters et al. (1980) 
Browne et al. (1981) 
INP-1 
INP-2 
INP-3 
INP-4 



20 



40 60 

E (keV) 



80 



100 



FIG. 6: Comparison between calculated and measured (open 
diamonds from [jj and open circles from fl4|) neutron capture 



cross sections on 



<Os. 



This leads us to reject the calculation 
187 Re(7,n) 186 Re reaction. 



■INP-1' for the 



In the particular case of the stable 187 Os target, direct 
experimental data are available for the 187 0s(n, 7) 188 0s 
reaction 0, and used as additional constraints on the 
nuclear ingredients, namely the combination of the nu- 
clear level density and optical potential. The resulting 
cross sections obtained with the 4 sets of nuclear inputs 
given in TablelTTIare shown in Fig^and seen to agree rel- 
atively well with experimental data, except in the specific 
'INP-3' calculation where the use of the global optical 
13 potential |39j | give rise to a cross section with an energy 
dependence relatively different from the one measured. 

The stellar Maxwellian-averaged neutron capture cross 
sections for the three reactions studied here are shown in 
Fig. [7| It is evident that although the photoneutron cross 
section is relatively insensitive to the nuclear level den- 
sity and neutron-nucleus optical potential, particularly in 
the 186 W(7, n) 185 W case, these quantities can affect the 
reverse rate significantly. If we characterize the remain- 
ing uncertainty affecting the prediction of the neutron 
capture rate by the ratio between the upper and lower 
limits obtained in Fig. [7| we find at an energy of 25 keV 
a factor of 1.9 for the 185 W(n, 7 ) 186 W reaction, of 2.1 
for 186 Re(?i,7) 187 Re and only 1.1 for 187 0s(n, 7 ) 188 0s. 
In the last case, the small error bars arise from the ad- 
ditional constraints made available through the experi- 
mental 187 0s(n, 7) 188 0s cross section. 
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TABLE III: Comparison between the F a values obtained in 
the present work (including INP-3 calculation) and those of 



Ref. 


|7j for different temperatures kT (expressed 


in keV). 


kT 


Present 


Ref. [7] 


12 


0.901 - 0.937 


0.867 


20 


0.879 - 0.886 


0.839 


25 


0.866 - 0.874 


0.830 


30 


0.859 - 0.867 


0.820 


52 


0.822 - 0.847 


0.813 



present calculation agrees relatively well with the value 
of 0.80 < F a < 0.83 at kT=30 keV from Ref. E3, but 
not with the value of 1 < F a < 1.15 from Ref. (10|. 

We now consider implications of the F a values con- 
strained by the present study to the Re-Os cosmochronol- 
ogy. The most tantalizing aspect of the Re-Os chronology 
is that it requires a rather detailed model of the chemical 
evolution of the Galaxy. Since constructing a reasonable 
model of the Galactic chemical evolution is beyond the 
scope of the present study, we rather focus on the uncer- 
tainty in the estimate of the age of the Galaxy within a 
schematic model. 

We recall here that abundances of elements in the 
relevant mass region can be symbolically expressed as 
follows: 



186 
187 
187 
187 



Os® 


= QgS _ 


h 186 Os p 


Os Q 


= 187 Os s - 


h 187 Os c 


Re Q 


= 187 Re r H 


- 187 Re s - 


Os c 


= 187 Re c 





7 Re c 



FIG. 7: Stellar Maxwellian-averaged neutron capture cross 
sections for 185 W (upper panel), 186 Re (middle panel), and 
187 Os (lower panel) calculated on the basis of the HF input 
defined in Table El 



V. IMPLICATIONS TO THE RE-OS 
CHRONOMETRY 

At stellar temperatures relevant to the s-process nucle- 
osynthesis (T ~l-3xl0 8 K), the 187 0s first excited state 
at 9.75keV is strongly populated and can significantly af- 
fect the estimate of the stellar neutron capture rate on 
187 0s. The correction to the cross sections due to the 
neutron capture on the 9.75 keV state is introduced by 
F(j in Eq. ©. 

On the basis of the present photoneutron data and 
the calculations, the temperature dependence of the F G 
factor has been re-estimated (Fig. [S]and Table ITTT|> . At 
the s-process temperature of 3 x 10 s K, all calculations 
converge to the F a value of about 0.87. However, at 
lower temperatures, the model INP-3 based on the global 
phenomenological optical potential predicts significantly 
larger F a values. This difference mainly originates from 
deviations seen in the laboratory cross section (cf. Fig.^jJ 
and should therefore be given a lower credibility. The 



Here, 0, s, p, r, and c represent the solar, s-process, 
p-process, r-process, and cosmoradiogenic origins, re- 
spectively. We introduce the following approximations: 



'Os Q 



186 0s p = px 

187 Re s = 

i 87 0s s = [F CT (a)( 186 0s)/( ( 7}( 187 0s)] 
approximation: Eq. (|T|)) 



186 



Os s (local 



Following the recent p-process calculations |44|, we es- 
timated p as lying within the 0.01-0.04 range, being in 
agreement with p=0.02 As noted earlier, we ignore 
the possible s-process contribution to 187 Re through the 
branchings at 185 W and 186 Re. 

Thus, the cosmoradiogenic component of 187 0s can be 
obtained as 



7 Os c = 



'Os 



(1-p) [F a (a) ( 18fa Os)/ (a) ( 187 0s)] 18b Ot 

(5) 

A consideration of the simplest model of a closed sys- 
tem would lead us to the assumption that the evolution 
of 187 Re can be effectively described by 



1187 



Re(t) 



dt 



-X e / fl87 Re(t) + Y{t) 



(6) 
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where A^ is the effective /3-decay rate of 187 Re in con- 
sideration of some enhancement by astration from the 
laboratory decay rate of Xp = ln2/(41.6 Gyr), whereas 
Y(t ) term represents the net r-process yield. As in Ref. 
|l9| . we adopt a simple form of Y(t) — y exp(— Xt) where 
A is a free parameter. Thus, we have 



187 



Re(t) 



Vie 



-Xt] 



A- A 



eff 



(7) 



Using d 187 Os c (t)/dt = + X 



eff 187 





Re(t), the abundance 
ratio between 187 Os c (t) and 187 Re (t) can be obtained as 



187 



Os c (t) 



187 



Re(t) 



B 
A 



(8) 



where 



1.00 
0.98 
0.96 
0.94 
[i^ 0.92 
0.90 
0.88 
0.86 





IMP H 




INP-2 




INP-3 




INP-4 : 



















1 2 3 4 5 

T [ 10 8 K ] 

FIG. 8: Prediction of the F a factor on the basis of the different 
HF calculations defined in Table UTI and compared in Figs. El 
and |S| with experimental data. 



A = e- x e ft - e~ xt (9) 



B = [1 - e _A 5 *] - [1 - e~ xt }X e / f /A. (10) 

The abundance ratio in Eq. JSJ at 4.55 Gyr ago is 
matched with that from Eq. J5J . 

For the metcoritic abundance of cosmoradiogenic 
187 Os c relative to 187 Re®, we used the following solar 
abundances pUf: 

187 Re &/186Q s Q = 3 51 ± o 09 

187 0s Q / 186 0s Q = 0.793 ± 0.001. 

The (n, 7) cross section ratios at s-process temperatures, 
cr( 186 0s)/er( 187 0s), were taken from Ref. Thus, 
the meteoritic quantity is determined with the present 
F a value being a unique parameter. As for X^ /Xp is 
used to calculate the abundance ratio in Eq. JSJ at t 
= Tq (the age of the Galaxy) -4.55 Gyr, Clayton H| 
derived 1.4 from the work of Yokoi et al. Q, whereas 
more recent analyses suggest considerably lower values 
|47) . We adopt here 1.2 as our standard value for the net 
enhancement by astration of the 187 Re /3-decay rate. 

Matching conditions in 187 Os c / 187 Re® were investi- 
gated in the Tq range of 11 - 15 Gyr, showing good 
agreement with < A < 2 Gyr -1 . We summarize some 
conclusions obtained under the simplest assumption of 
r-process nucleosynthesis yields varying exponentially in 
time. 

• F a values at typical s-process temperatures are in 
the range of 0.86-0.94 (Table |TTTJ|. 

• The probable range of the differential coefficient 
dT G /dF a is -(5.0-12.8) Gyr. 



• Consequently, the remaining uncertainty of Tq 
that stems from that of F a values is less than 
1 Gyr. When the temperature dependence of 
(cr)( 186 0s)/(cr)( 187 0s) is considered along with 
that of F a , the uncertainty in Tq is approximately 
halved. 

We note here that the model of Yokoi et al. [jjj can- 
not be reconciled with the present F a data. The model 
of chemical evolution developed there favored F a values 
much higher than unity. If we mimic the results in terms 
of Eq. the corresponding values of A become nega- 
tive. This also explains the much larger \dTQ/dF a \ values 
of up to 100 Gyr (as seen in the slope of Fig. [5] in Ref. 
0). Finally, it' is noted that if we would use the age of 
the universe as derived by the WMAP (T v = 13.7 ± 0.2 
Gyr) 0|, and assume T v - T G m 1 Gyr, A values that 
are consistent with the present F a value are narrowed in 
the ranges of 0.16-0.46 Gyr" 1 , inclusively, and 0.22-0.34 
Gyr -1 , exclusively. The exclusive values are commonly 
allowed in calculations with all possible combinations of 
kT, F a , and Tq. 

VI. CONCLUSION 

Photoneutron cross sections were measured with 
accuracy for 186 W, 187 Re and 188 Os using quasi- 
monochromatic photon beams from laser Compton scat- 
tering (LCS) at energies near the neutron thresholds. 
The cross sections were used to constrain the model 
parameters in the framework of the Hauser-Feshbach 
model. Four different sets of nuclear ingredients were 
adopted to estimate the photoneutron cross sections 
and the reverse radiative neutron capture cross sections. 
When no experimental data on the direct (n, 7) cross sec- 
tion is available, an accuracy of about a factor of 2 was 
achieved in the predictions. The influence of the neutron 
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capture by the 9.75 keV first excited state in 187 Os which 
is substantially populated in stellar plasmas at typical s- 
process temperatures has been estimated in connection 
with the 187 Re- 187 0s cosmochronology and shown to lead 
to an increase of the neutron capture rate by a factor of 
about 1.15 at a temperature of 3 x 10 8 K. Uncertainties 
by about 10% associated with the neutron-nucleus opti- 
cal potential still affect the stellar rate at temperatures 
between 1 and 2 x 10 8 K. 

The correction factor F a to be used in the local ap- 
proximation (Eqs. and @) was constrained well in 
the present study fTablc ITTT|) . Based on the simplest as- 
sumption of r-process nucleosynthesis yields varying ex- 
ponentially in time, the cosmochronological uncertainty 
in the age of the Galaxy arising from the F a values is 



estimated to be less than 1 Gyr; when the temperature 
dependences of both (cr}( 186 0s)/(er)( 187 0s) and F a are 
considered, the uncertainty is less than 0.5 Gyr. 
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